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Abstract: We investigate the enhancement of neutraUno rehc density in the context of a 
reahstic cosmological scenario with Quintessence. The accurate relic density computation 
we perform allows us to be sensitive to both cases with shifts in the abundance at the 
percent level, and to enhancements as large as 10^. We thoroughly analyze the dependence 
on the supersymmetric spectrum and on the mass and composition of the lightest neu- 
tralino. We point out that supersymmetric models yielding a wino or higgsino-like lightest 
neutralino become cosmologically appealing in the presence of Quintessence. 
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1. Introduction 

The latest years have been marked by fast and steady progresses in observational Cosmol- 
ogy. The picture emerging from the data is remarkably self-consistent [||, and, at the same 
time, it is significantly different from that which used to be the standard lore less than a 
decade ago. The greatest surprise came with the discovery that the Universe is accelerat- 
ing [^, with most of its energy density currently in the form of "dark energy" , a component 
with negative pressure. In the simplest scenario compatible with these observations, with a 
cosmological constant term which has very recently (on a cosmological timescale) matched 
and overtaken the rapidly decreasing matter term, a large degree of fine-tuning is needed in 
the initial conditions in the early Universe. In recent years, alternative models have been 
proposed, in which dark energy is introduced as "Quintessence" |^, ^, P], a dynamical 
term, such as a scalar field rolling down an effective potential, gravitationally coupled to 
the other components of the Universe and possibly tracking the largest of them long before 
becoming the dominant one at the present epoch. Besides driving the dynamics of the 
Universe today, such a term, which has a time varying equation of state, may have played 
a role during other stages in the evolution of the Universe: other issues in the "canonical" 
picture, not directly tested by observations so far, might have been overlooked, a problem 
which has been addressed in a number of recent analysis and in different contexts, see, 
e.g., Si. 

The latest Cosmological observations have confirmed that there is a large discrepancy 
between the total amount of matter in the Universe and the baryonic component; in terms 
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of the ratio 0, between the mean density and the critical density, the contribution of non- 
relativistic matter has been found Q to he in the range Jlm^^ = O-l^Slggog (here h is the 
Hubble constant in units of 100 km s~^ Mpc^^; h = O.Tl^QQg [|l|), much larger than that of 
the baryonic component i^i^h'^ = 0.0224 it 0.0009. The nature of cold dark matter remains 
one of the main puzzles in today's research. Probably, the most natural approach to this 
problem is to suppose that dark matter is just another of the thermal leftovers from the 
Early Universe. It is a remarkable coincidence that the thermal relic abundance of stable 
weakly- interacting massive particles (WIMPs) is roughly of the order of Qm- in several 
classes of extensions to the standard model (SM) of particle physics such dark matter 
candidates arise naturally, the lightest neutralino in supersymmetric extensions to the SM 
being the most widely studied case. 

As in any process involving departure from thermal equilibrium, the freeze-out of 
WIMPs can be understood in terms of a particle interaction rate F falling below the 
expansion rate of the Universe H. T and its scaling with temperature are fully defined 
once the particle physics setup has been specified; what, instead, about H and its thermal 
scaling? The freeze-out of WIMPs is predicted to take place at a temperature Tf.^ in 
the range of a few GeV or above, much earlier, in the history of the Universe, than any 
of the processes which can be directly tested through cosmological observations. The 
standard prejudice is to extrapolate the snapshot of the Universe we derive from tests of the 
standard big bang nucleosynthesis (BBN) , at a temperature of about 1 MeV, backward to 
the WIMP freeze-out temperature, assuming a radiation dominated Universe. Implications 
of nonstandard cosmological scenarios at the time of decoupling were first envisaged and 
analyzed by Barrow and by Kamionkowski and Turner |l^: the various possibilities 
ranged from setups where a shear energy density component is generated by anisotropic 
expansion |9|, 10 1 to Brans-Dicke- Jordan cosmological models [10|. Following that line, in 
a recent analysis, Salati pointed out that predictions for the WIMP relic abundance 
are significantly different if, on the other hand, one supposes that the freeze-out happens 
during a period of Quintessence domination, as it can be the case if the Quintessence field 
undergoes a phase of "kination" (see the next Section for a clarification on this point). 

The underlining idea is rather simple: if on top of the radiation energy density a 
Quintessence contribution is added, the Universe is forced into having a faster expansion. 
A larger H implies that the matching between annihilation rate and expansion rate takes 
place at a higher temperature, hence the WIMP equilibrium number density freezes in 
at a larger value, giving a net increase in the final WIMP relic abundance (even a mod- 
est increase in Tj induces large effects, because the freeze-out of WIMPs happens when 
WIMPs are non-relativistic and their equilibrium number density is along the exponentially 
suppressed Maxwell-Boltzmann tail) . We examine here this mechanism in detail: we imple- 
ment a realistic prototype for Quintessence into the density evolution equation describing 
the decoupling of a WIMP dark matter species, and solve such equation numerically with 
the state of the art technique as developed in the DarkSUSY package [0, |l^, |l^ . We are 
then in the position of giving some general criteria to quantify the quintessential enhance- 
ment of relic densities, and of illustrating what classes of models are more sensitive to such 
an effect. We present results for a few sample cases of neutralino dark matter candidates 
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in the context of the minimal super symmetric extension of the SM (MSSM). In particular 
we discuss at length the case of higgsinos and winos, which are usually disregarded as dark 
matter candidates, as their thermal relic abundance is typically exceedingly small in the 
standard cosmological scenario. 

The outline of the paper is as follows: In the next section we introduce our setup for 
quintessential Cosmology, and in section 3 we illustrate the possible role of Quintessence 
in shifting the neutralino relic abundance. In section 4 we focus on the case of higgsinos 
and winos, while section 5 concludes. 



2. A role for Quintessence in the early Universe 

We consider a Friedmann-Robertson- Walker cosmology including, in addition to the usual 
radiation and matter terms, a Quintessence component, which we sketch as a single spatially- 
homogeneous scalar field (j) with a potential V {(p). The energy density and pressure asso- 
ciated to this component are, respectively: 

1 /d(hV „ , 1 



By formally writing the equation of state for (j), = wptf,, we see that w can vary between 
and —1, going from the regime in which the kinetic term is much larger than the 



potential term ( "kination" phase [ 14 , |6[ ) to the opposite case when the field is frozen into 
one configuration and just behaves as a cosmological constant term. 

To trace the Hubble expansion rate H as a function of the temperature T, we need to 
solve the Friedmann equation (we assume a spatially fiat Universe): 



1 

3M^i 



(T) = 1^ [Pr (T) + Pm (T) + (T)] , (2.2) 



coupled to the equation of motion for 

d^(j) d(h dV 



In Eq. (O) above, Mpi = l/VSvrG is the reduced Plank mass, while pr and pm are. 



respectively, the energy density in radiation and in matter. To find Pm{T) we simply scale 
its current value with the appropriate equation of state; for Pr{T) the analogous procedure, 
often implemented in Quintessence studies, is not accurate enough for our purpose. We 
write instead Pr{T) as: 

Pr{T) = ^gMT)T' (2.4) 

and compute g^siT), the effective degrees of freedom coefficient, summing the contributions 
from all particles in the context at hand (i.e. , in the super symmetric extension to the SM 
we consider below, we find ges(T) for each of the mass spectrum we generate, including 
both standard model particles and the supersymmetric partners; details on how to compute 
9es{T) are given, e.g., in Ref. |5|). 
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P(j,{T) and H{T) can be derived once we specify the initial conditions and the potential 
V {(j)). The set of initial conditions we resort to is not the most generic in a context of 
quintessential cosmology, as we wish to restrict to solutions in which is initially larger 
than pr, and then is red-shifted away more rapidly than the radiation component (this 
happens in the kination phase) so that the radiation domination epoch can start before 
the time of BBN. This corresponds, e.g., to the reheating scenario suggested in Ref. |16| 
and, more generically, in models that try to generate, through the same mechanism, the 
current inflationary period and a period of inflation in the early Universe (see, e.g., the 
quintessential inflationary models of Ref. [0]). As regards the potential, we choose to work 
with the exponential form Q: 



V {(f)) = Ml, exp (-A(/./Mp) 



(2.5) 



one of the simplest examples of potential which can lead to an attractor solution, self-tuning 
the contribution of 4> to the energy density to the "background" contributions, i. e. , in our 
case, radiation plus matter. We will select regions in the parameter space in which the 
attractor is reached, so that we deal with a Quintessence model playing a main dynamical 
role in the recent past; actually the simple form for the potential we picked is not going to 
give the behavior of observed today (the solution converges to the same equation of state 
as the largest background component), however slight modiflcations to it can introduce the 
right scaling, see, e.g., Ref. |18, On the other hand, we checked in a few test cases that 
this slight reshuffling of the potential can be introduced without varying signiflcantly the 
dynamics of the field in the early Universe^; we can therefore safely ignore details in the 
late time behavior of the field. The exponential form itself is not critical for the results 
we will present, however, in order to have a viable prototype for Quintessence, one should 
restrict to setups which allow to reach tracking after going through a kination phase, a 
feature which is not shared by all potentials which do admit an attractor solution; further 
examples of such setups are given in Ref. ||2^ . 

The possibility of having a kination regime with the exponential potential introduced 
above was examined in detail by Ferreira and Joyce in Ref. Q; we closely follow their 
analysis here, choosing, in particular, the same set of initial conditions. In terms of Hi, 
the expansion rate at the end of inflation, we flx: 



6init) = 3Mj,Hf , ^ 



, and pr (Ti 



init } 



(2.6) 



the energy scale for Hi can be in the range [10*^,10^^] GeV Q^. For any given choice of 
the two parameters in our model, i.e. the parameter A in the exponential potential and 
Hi, the system of differential equations ( |2.2| ) and ( |2.3D can be easily solved numerically. 
A few sample cases (see also Ref. |^]) are shown in Fig |l] where we plot energy density 

^We explicitly worked out the cases of the hyperbolic cosine and of the so-called AS potential, i.e. 
an exponential potential with a power-low prefactor |l8j. 

^From WMAP data a more stringent, though model dependent, upper bound has been derived 
Hi < 3 X lO^^GeV, possibly cutting out a slice in the parameter space. This, however, does not affect the 
regime of interest in the present paper, which lies in the low Hi range. 
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Figure 1: Sample cases of Cosmologies with a Quintessence field in a kination phase dominating 
the total energy density in the early Universe. The solid lines refer to the sum of the radiation 
and matter energy density components, the other curves to the Quintessence contributions. We 
have considered the case of an exponential potential and show here the effect of varying the two 
parameters in our setup: In the left panel three sample values for the expansion rate at the end of 
inflation Hi are chosen, in the right panel the parameter A in the exponential potential is varied. 

components as a function of the Universe scale factor a. Solid lines refer to the sum of 
the matter and radiation component (with an ankle at the value when the former, scaling 
as a~^, takes over the latter, which scales instead as a~'^); Quintessence energy density 
components (scaling as in the initial kination phase) are also shown. An increase in 
Hi shifts the initial conditions to a higher temperature and, at the same time, lowers the 
initial ratio p^j Pr ^ ^T^'- the effect is obviously to anticipate the radiation dominated 
phase and the tracking phase. An increase in A makes the potential steeper, henceforth 
the kination phase to last longer; from the point of view of the early time behavior, there 
is just a slight decrease of the normalization of in the kination phase due to a quicker 
transition from the initial condition (with (j) at rest) into the fast rolling phase (for this 
specific potential the normalization of p^ at tracking scales instead as 1/A^ Q, but, as 
already mentioned, we can neglect here details in the late time behavior of p,^). 

In Fig I we show the region of the parameter space Hi versus A which is compatible with 
the limits from BBN. The latter are often reported in terms of bounds on extra relativistic 
degrees of freedom AA'eff at BBN, which may be translated in an upper bound on the 
fraction il,^ = p^p/pr ^ 0.1 for AA'^eff — 1-0. In the dark shaded region the Quintessence 
field produces a contribution to the energy density at 1 MeV which exceeds the stated 
constraint. In the lower part of the figure this happens because the (j) field is still in the 
kination phase, and exceeds at the time of BBN the bound < 0.1; in the upper- left 
part, instead, the field is already in the attractor solution at 1 MeV, and again violates 
the BBN bound. The light shaded region indicates, on the other hand, the cases in which 
the ratio between p^ and the background components gets overly damped, so that tracking 
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Figure 2: Region of tlie parameter space Hi versus A wliich is excluded by the limits from standard 
big bang nucleosynthesis. Also indicated is the region where the ratio Quintessence energy density 
to background energy density gets exceedingly small and tracking does not take place before the 
present epoch. 

does not take place before the present epoch {overshooting solutions); both this cases are 
disregarded in the present work. 

3. Quintessence and neutralino dark matter 



The density evolution equation describing the decouphng of a species from the thermal 
bath is the Boltzmann equation: 

^ = -3i?n - (cTcffv) (n^ - n^q) , (3.1) 

where n is the particle number density, neq the corresponding equilibrium density, and 
(cjeffi') the thermally-averaged pair annihilation cross section. In case coannihilation effects 
p^ , p^ , p^ ] are present, i.e. when extra particles are nearly degenerate in mass with the 
(lightest) stable one, and share a quantum number with it, n is defined as the total number 
density summed over all coannihilating particles, and {(Jesv) contains a weighed sum over 



all relevant rates, see, e.g., ||2J] for details. Rephrasing the density evolution equation (^J) 
in terms of y = n/s, where s is the entropy density, as a function of the ratio between the 
particle mass and the temperature, x = m/T, one gets the equation 



dY 

dx 



3H dx 



This equation can be recast as 
X dY 



Veq dx 



T X ds 
H 2>s dx 



y2 
^ eq 



1 ~ - — 



r 

H 



y2 
^ eq 



(3.2) 



(3.3) 
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where the roles of the annihilation rate T = neqicr^sv) and of the expansion rate H become 
more transparent: When T/H ^> 1, thermal equilibrium holds, that is y ~ Y^q, while in 
the opposite regime Y^q <^ Y. As a rule-of-thumb, the freeze-out occurs when the ratio 
T/H ~ 1. An increase in {(icfiv) yields a decrease in the freeze-out temperature. On the 
other hand, as already mentioned, if we increase H adding an extra component on top of 
the radiation component, such as in the quintessential scenario outlined in the previous 
section, the freeze-out takes place at larger temperatures, and the thermal relic density can 
be sharply enhanced. Even a modest shift in Tf.^ has a dramatic effect, with an increase 
in the relic density of typical dark matter candidates which can be up to several orders 
of magnitude: This is due to the fact that the thermal equilibrium number density scales 
exponentially with the temperature, i.e. Ugq oc (mT)^^^ exp(— m/T). 

To better quantify the enhancement, rather than introducing a simplified toy- model 
and deriving an approximate solution to Boltzmann equation as done in Ref. 0, we in- 
terface the expansion rate of the Universe H we have derived in the Quintessence scheme 
described in sec. 2 into the accurate numerical solution of Boltzmann equation included in 
the DarkSUSY software package [12, 11, 13 1. For any given set of initial conditions in the 
Early Universe, and for any particle physics setup with given mass spectrum and particle 
couplings, we then have a tool to compute WIMP relic abundances with an accuracy of 
the order of 1%, i.e. at the precision level of upcoming measurements of f^m- 



3.1 The particle physics setup 

In our numerical study we focus on neutralino dark matter candidates within the context 
of the MSSM. We recall that neutralinos are mass eigenstates derived from the linear 
superposition of neutral gauginos, the bino B and the wino W^, and higgsinos: 

Xi = NaB + Ni2W'^ + iV,3^i° + N,^H^ . (3.4) 

The lightest of them is often also the lightest supersymmetric particle (LSP), and, thanks 
to i?-parity, stable; such neutral and stable particle is potentially a good WIMP dark 
matter candidate. We define wino fraction of the lightest neutralino the quantity |A'^i2| and 
higgsino fraction the quantity y^jiVisP + |-/Vi4p. If the lightest neutralino is wino-like, i.e. 
if it has a dominant wino fraction, it is also approximately degenerate in mass with the 
lightest chargino; analogously, if it is higgsino-like, there will be a quasi-degeneracy both 
with the next-to-lightest neutralino and with the lightest chargino. These degeneracies 
yield large coannihilation effects. Particular models which entail a wino or higgsino-like 
LSP will be discussed in sec. 

As shown in eq. ( |3.2| ), the neutralino relic density strongly depends on the neutralino- 
neutralino annihilation rate. It is well known that, in the case of bino-like LSPs, this 
rate tends, in general, to be not large enough to avoid an overproduction of thermal relic 
neutralinos. It is often quoted that, within the so-called constrained MSSM, only four 
regions of the parameter space give rise to relic densities compatible with the current 
cosmological upper bound on Qrn- (1) A bulk region, with a light neutralino and, typically, 
tight accelerator constraints {e.g., BR{b — > 57), (5a^ or the bound on the mass of the 
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Figure 3: Allowed region and iso-Hi curves in the (A, Afi) plane. The dark blue region does not 
fulfill BBN constraints, while in the light blue zone no tracking takes place before nowadays cosmo- 
logical times (overshooting) . We show with bold lines the curves with the same =, 10^, 10"'^'', 10^^ 
and 10^^. Also shown, in dotted lines, curves at equal logarithmic spacing from those in bold lines, 
i.e. at around 1.58, 2.51, 3.98 and 6.31 times the upper bold line Hi. 



lightest CP-even Higgs boson); (2) A coannihilation region with the NLSP, which in most 
of the parameter space is the lightest stau; (3) A funnel region, where the neutralino mass 
is close to half the mass of a neutral Higgs boson, most often the CP-odd Higgs boson 
A, entailing strong resonance effects which drastically enhance the neutrahno-neutralino 
cross section; (4) A focus point region, at large mo, where electroweak symmetry breaking 
conditions dictate a sizeable higgsino fraction in the LSP composition. 

In a more general setting, when one relaxes some assumptions of the constrained 
MSSM in the context of supergravity, or resorts to other phenomenological supersymmetry 
breaking mechanisms, the lightest neutralino may naturally feature large higgsino or wino 
components. For instance, if the gauge kinetic function depends on a non-singlet chiral 
superfield, whose auxiliary F component is responsible for the gaugino mass generation, 



particular (fixed) ratios exist between the three gaugino masses at the GUT scale |2^, 26 1. 
In the simplest case of SU{5) GUTs, the three non-singlet representations are the 24, the 
75 and the 200. The latter two (largest) representations yield higgsino-like LSPs. In what 
follows, we will pick as a sample case the 200, which has the virtue of allowing a larger 
parameter space, particularly at large tan/?, with respect to the 75. 

3.2 Relic abundance enhancement 

We choose as indicator of the enhancement produced by a quintessential cosmology onto 
the neutralino relic density the quantity 
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Figure 4: Ratio Quintessence energy density to radiation energy density as a function of temper- 
ature, for a few sample choices of the paramters A and Hi. 



O/i^ being the neutralino reUc density in a cosmology without Quintessence (standard case), 
and ^qh"^ that with Quintessence. In Fig ^ we show a first example of how large an 
enhancement we can obtain in our Quintessence setup. In the (A, Ail.) plane, we plot curves 
at fixed Hubble rate at the end of inflation Hi ; in this figure we have focussed on a specific 
neutralino candidate, selecting the point in the bulk region of minimal supergravity defined 
by the supersymmetric parameters M1/2 = 2300 GeV, tuq = 3000 GeV, tan (3 = 45, Aq = 
and fi > 0. With such choice we avoid, for the moment, coannihilation or resonance effects, 
and consider the case of a heavy (bino-like) neutralino {rriy, ~ 1 TeV), for which the 
enhancement effects are maximized: as we will discuss in detail below, the heavier the 
freezing-out particle, the larger the quintessential effects at decoupling, and therefore the 
larger Afl. The shadings in the figure reproduce those of Fig dark blue corresponds 
to a scenario ruled out by BBN constraints, while light blue indicates overshooting (i.e. 
tracking is not achieved at current cosmological times). The largest enhancement allowed 
is of the order 10^ and it is about the largest enhancement we can get with a WIMP of 
mass smaller than 1 TeV. Moreover, we see that the minimum variation we are sensitive 
to, an enhancement of the order of 1% occurs for Hi slightly above 10^^: This is the limit 
around which, in the process of decoupling of a WIMP with mass below 1 TeV, we are not 
sensitive any longer to the presence of a kination phase. 

As a last remark on Fig|^, we comment on the shape of the iso-Hi curves at low A: the 
abrupt raise in AQ is due to the fact that tracking takes place, at low A, as soon as /j,^ ~ p^j 
without an overshooting period where p,^ <C Pr- This of course yields a large quintessential 
contribution to H at all time scales relevant to neutralino freeze-out, contrary to the case 
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where an overshooting {p^ <C Pr) takes place. However, these sudden rises in Afi happen in 
the region excluded by BBN constraints, as the tracking phase cannot appear before BBN. 
Rejecting these cases, we find, as expected, that all large enhacements take place for models 
in which the field i;^ is in a kination phase, p^ oc . It appears then natural, rather than 
searching for a parametrization of Ail as a function of A and to introduce, as one single 
relevant parameter, the relative normalization of pcf, and pr at a given temperature. We do 
not have the freedom to choose such temperature as low as, say, the BBN scale T ~ 1 MeV, 
as it is not unusual for the field (j) to be, at this temperature, in its transient phase between 
kination and tracking (this is true even for our Quintessence model in which the transient 
is always very sharp, with (j) behaving like a cosmological constant; more generic models 
allow for softer and more complicated transients). This is implicit in Fig |^ and more 
directly illustrated in Fig ^ where the ratio Pcj,/ Pr is plotted versus the temperature for 
a few sample choices of the parameters A and Hi: in the kination phase this ratio scales 
nearly as (behaviour of all curves for the largest temperatures display), but this is not 
necesseraly true at 1 MeV. 

In order to quantify the amount of Quintessence relevant for the neutralino relic density 
enhancement, we choose to resort to the parameter: 



- ^ (T/:^j , (3.6) 

where we defined the temperature of neutralino freeze-out in absence of Quintessence 
(we use here the convention that the freeze-out temperature is the temperature at which 
the abundance of the relic species is 50% larger than the equilibrium value, as computed 
from the full numerical solution of the density evolution equation) . Avoiding the reference 
to a single cosmological temperature, the parameter entails a nearly model-independent 
estimate of the relevant amount of Quintessence. In order to have a rough estimate of this 
effect, we can consider as a first order approximation of Vt the following expression: 

Pc Pc ~ Pc s{Tf,,,) Pc s(Tf.„.) (acsv) 

where in the third step we introduced the approximation that Y(To) ~ y^q(Tf ^ ), while 
in the last step we assumed that the freeze-out temperature is defined by the relation 
T{T,,J = H{T,,J. 

In the context of a quintessential cosmology, the expansion rate gets an additional 
factor given by: 



\ Pr V (-^f.o. j 

where the last approximate equality holds under the simplifying hypothesis that p^ oc a~^. 
Since the actual shift in the freeze-out temperature should not be dramatic, i.e. T^^ ~ T^^, 
we expect the relic abundance enhancement to be roughly driven by the factor y^l + C^- 
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Figure 5: The relic density enhancement Aft as a function of the Quintessence-to-radiation ratio 
at neutrahno freeze-out for four SUSY models with neutralino masses equal to 200, 500, 750 and 
1000 GeV. In all cases, tan/3 = 45, /i > and Aq = 0. The parameters mg and M1/2 are tuned in 
order to obtain the desired neutralino mass and a mass splitting Attia = {2m^ — m a) / {'i-'niy^) = 45%. 



3.3 Wandering about the SUSY parameter space 

As already alluded, the phenomenon of quintessential enhancement of the neutralino relic 
density should depend not only on the particular quintessential setup, but also on the 
details of the choice of the supersymmetric dark matter candidate. In particular, we 
will show here that there is a dependence (1) on the lightest neutralino mass, since this 
sets the overall energy scale at which neutralinos freeze-out, (2) on the details of the 
composition of the lightest neutralino in terms of its bino, wino and higgsino fractions, and 
(3) on the supersymmetric mass spectrum, especially if it entails coannihilation processes 
or resonances. 

We focus first on the role of the neutralino mass. In the left panel of Fig ^ we plot 
Ail versus the ratio of the Quintessence energy density to the radiation energy den- 
sity computed at the freeze-out temperature, for four neutralinos in the bulk region of 
the minimal supergravity model and with masses nriy- = 200, 500, 750 and 1000 GeV. 
The supersymmetric parameters have been chosen in order to exclude coannihilations or 
resonance effects, and the neutralinos are always almost completely bino-like. As regards 
the cosmological setup, we fixed A = 3.5, and let Hi vary in order to span over S^^. As 
already mentioned, larger masses yield larger enhancements for a given Quintessence setup; 
we can deduce however from the figure that, to a good approximation, for relatively low 
quintessential fractions, the quantity fixes Afi. In particular we find that a 10% en- 
hancement in the relic abundance is obtained if at the Quintessence energy density is 
about 60% of the radiation energy density, and that a 100% effect is generated by ~ 12. 
As the amount of Quintessence grows, the spread in enhancements for different masses, 
and at a fixed ^(^, starts to increase, as shown in the right panel of Fig ^; AO is in the 
range 60 to 80 for the largest value of the ratio we display, = 10^, a 15% spread around 
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Figure 6: The quintessential enhancement of the neutraUno rcHc density for the four benchmark 
models of Tab|| as a function of respectively at large in a linear scale (left) and at small 
in a logarithmic scale (right). 



the mean value: this reminds us that the freeze-out is not a sudden process characterized 
by a single energy scale, and if an accurate estimate of relic abundance is needed, a careful 
numerical solution to the Boltzmann equation has to be implemented. 

In order to study the dependence of the quintessential enhancement on the particular 
structure of the lightest neutralino, we focus on four benchmark points, chosen ad hoc to 
have the same freeze-out temperature in absence of Quintessence — 9.35 GeV. We 
report in Tab. |l| the details of the four benchmark models under scrutiny. The model 
with the label bulk is characterized by a mass splitting with the A Higgs boson AniA = 
(2m^ — ma)/{2m^) ~ 45% and with the NLSP, which is the lightest stau, of Am^^ = 
(mf^ — m^/m-^ ~ 85%, and is therefore outside the resonance or coannihilation regions. 
The model in the stau coannihilation region achieves nearly complete mass degeneracy 
between the neutralino and the lightest stau, Am-f^ ~ 0, while the model in the funnel 
region has Am^ ~ 0. Finally, the model representative of the higgsino dark matter case is 
picked in the representation 200 and has a higgsino fraction larger than 95%. 

We analyze in Fig |6| the quintessential relic density enhancement generated in the 
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0.009 
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Table 1: The four Benchmark models representing different possible scenarios for the lightest neu- 
tralino. The scalar trilinear coupling = in all cases, and the neutralino freeze-out temperature 
in absence of Quintessence is, again for all four scenarios, fixed at rf„^ ~ 9.35 GeV. 
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Figure 7: {Left ): The effective anniliilation cross section Wcs / {'iE'^g) as a function of poff for tlic 
four benchmark models of Tab. 1. The dashed line indicates the thermal weight factor K(poff , T) at 
the freeze-out temperature T^^. (Right): The thermally averaged effective cross section (fToffw) as 
a function of the temperature for the four benchmark models. The brown dotted lines represent the 



suppression factor (^^1 + C</>^~ j ) due to a quintessential component whose amount at freeze-out 
is given by the corresponding reference scale is shown on the right-hand side of the figure. 



four benchmark cases as a function of the Quintessence-to-radiation ratio at the neutralino 
freeze-out temperature without Quintessence (^(^). We let vary from zero up to 10^, 
highlighting the large and small regimes, respectively, in the left and right panels. We 
find that the region of the parameter space where the enhancement is more effective is the 
coannihilation strip. Between the coannihilation and the bulk regions, we find the interme- 
diate case of the higgsino LSP: this is in fact expected, since coannihilations phenomena 
take place there as well, though on top of neutralino pair-annihilation processes which 
are by far more efficient than in the bino-like case. Finally, the weakest quintessential 
enhancement is found in the funnel region. 

With the purpose of understanding the peculiar pattern emerging from Fig |6|, we 
consider the effective thermally averaged annihilation cross section (cTefff), which we rewrite 
in the following form p3l: 



(cTeSv) 



(3.9) 



where Peff and E'eff = \j + "^x respectively the three- momentum and the energy of 
a pair of lightest neutralinos in their center of mass frame. Weg is the effective annihilation 
rate per unit volume and unit time; in case coannihilation processes are present, it can be 



written as a weighed sum over all annihilation and coannihilation processes |24|: 



9\ 



(3.10) 
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Figure 8: A comparison between the thermally averaged cross section (acsv) as a function of 
the temperature, compared with the temperatures of freeze-out Tf „. (nearly vertical lines on the 
right) and of freeze-in Tf i_ (left lying lines). The first one is defined as the temperature when the 
neutralino number density becomes 50% larger than the equilibrium number density p^ , while 
Tf i is the temperature corresponding to the moment when the number density normalized to the 
entropy density gets 50% larger than its asymptotic value, i.e. Yq, which gives the nowadays relic 
density. 



Here gi is the number of internal degrees of freedom for the particle i, pij the common 
magnitude of the three-momentum in the center of mass frame for the process involving 
the particles i and j, aij the relative cross section, while the index 1 refers to the LSP (i.e. 
Wes correctly reduces to the annihilation rate Wn for a pair of lightest neutralinos, if no 
coannihilation is present). Finally, eq. ( p.9| ) defines the function K{pes,T), which contains 
the temperature dependence in the Maxwell-Boltzmann approximation. 

We plot in the left panel of Fig |^ the resulting Wes/{4:E'^g) for the four benchmark 
models of Tab. |^. When a coannihilation channel becomes kinematically allowed, the effec- 
tive cross section abruptly raises, as evident from both the case of the stau coannihilations 
and that of the higgsino LSP. In this last case, also notice how the cross section at low peS 
is larger than in the bulk, or stau coannihilations, case. The figure also shows the typical 
bell shape of the function k, which exponentially suppress the large Peff tails of Wcs upon 
integration over P(,g. As T is reduced, the peak of the corresponding k function is shifted 
towards lower pefT values. In the right panel of Fig |^ we plot the result of the integration 
over Peff) obtaining the corresponding curves as functions of T for the various cases. The 
shape of the resulting curves comes not as a surprise: a decreasing function of PeS (e.g. the 
case of the funnel) gives rise to a decreasing function of T, and vice-versa. 

The temperature evolution of ((Xcfr^) critically enters the thermal history of neutralinos, 
which, we remind the reader, proceeds from higher to lower temperatures, i.e. from right 
to left in Fig At very large temperatures, the number density of relic species closely 
follows the equilibrium number density. Afterwards two steps can be identified. First, 
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the number density n starts to b ecome largcT tlicin the eQuihbrium number density Ti^q^ 
in a stage commonly called freeze-out. This step takes place when the quantity F = 
n{aesv) becomes of the same order of the expansion rate of the universe H{T). (We 
recall that we did not define the freeze-out temperature Tf.o from the relation T = H, 
but rather as the temperature when the neutralino number density becomes 50% larger 
than the equilibrium number density, i.e. n = 1.5ncq Then, a second step goes on, 

where (co-) annihilations of neutralinos further reduce the number-to-entropy density Y 
towards its asymptotic value Yq, with an efficiency which depends on the relevant scattering 
cross sections at temperatures below freeze-out. For definiteness, we label as freeze-in 
temperature Tf the temperature at which Y = 1.5 Yq. 

In our discussion we have so far considered just the first step, i.e. the increase in Tf ^ 
set by the shift — > ff^Q 1 + P(j,/pr- To explain the pattern emerging from Fig 0, we 
need however to consider also the second effect mentioned above. The evolution of Y after 
the freeze-out temperature is described by the Boltzmann equation ( |3.2| ) with l^q = 0. 
Once exactly integrated, one obtains: 



Y„ r,.. Jr„ H 3dT 

In presence of Quintessence, the integrand in the right-hand side gets the suppression 
factor, which similarly to step 1 has the form 



111 11 



(3.12) 



Clearly, the longer the post-freeze-out annihilations last, the larger the suppression 
of the relic abundance is. The convolution between (aesv) and the suppression factor of 
eq. ( |3.12| ) sets the shift in the temperature range where neutralino annihilations are effec- 
tive, down to the freeze-in temperature. The four benchmark models show very different 
low temperatures behavior of (aesv): in the case of the funnel, {(Jesv) increases at lower 
temperatures, hence we expect a sharp shift in Tf ; to lower values as increases; on the 
other hand, in the case of stau coannihilations, {ctesv) rapidly drops to low values at low 
temperatures, therefore freeze-in takes place earlier, at higher temperatures. The higgsino 
LSP and the bulk cases are intermediate, and quite similar to each other. The emerging 
picture is clarified in Fig |8|, which summarizes the two effects we discussed: we depict both 
the freeze-out and the freeze- in temperatures as functions of which is shown on the 
vertical scale at the right-hand side of the plot, as well as the thermally averaged cross 
sections in the four cases under inspection. While the (slight) hierarchy in the freeze-out 
temperature depends on the cross sections at temperatures larger than freeze-out, the above 
mentioned convolution dictates the shape of the freeze-in curves. The pattern emerging 
from the superposition of the two mentioned effects is the one we found in Fig 

Since both the freeze-in and freeze-out temperatures depend on the same factor, which 



scales with as y 1 -\- £,^T'^ /T^^ , we expect to be able to give a rough rule-of-thumb 
estimate of the quintessential enhancement effects for a given We do not, however. 
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expect a rigorous scaling as the two intervening effects act on different temperature ranges. 
Tlie approximate scaling we propose tries to keep track of both the correct limit as — > 
and the would-be behavior at large ^(^: 

n^/n^ + 6-0.5 (3.13) 

We find that the actual enhancement in the full numerical computation reproduces, at large 
^cf,, our simple ansatz, although with a less steep increase: a best fit procedure gives in fact 
0-4 ^ b < 0.5. Moreover, the interplay between the shifts in the freeze-out and freeze-in 
temperatures (see Fig non-trivially dictates the value of the a parameter. From the 
overall fit to the benchmark model data we obtain 0.1 < a < 0.2. Though the accuracy 



of resorting to such an approximate formula is necessarily low, eq. ( 3.13 ) gives the correct 
order of magnitude of the quintessential enhancement, and may be of relevance in some 
cases. For instance, if one deals with a cold dark matter model giving rise to extremely low 
relic densities, and embeds it in a quintessential cosmological scenario, the formula given 
in eq. ( p.l3 ) may indicate the order-of- magnitude amount of Quintessence at freeze-out 



which the given model would need in order to get the correct relic density today. 



4. Resurrecting higgsino and wino dark matter 

When the lightest neutralino is dominated by its higgsino or wino component, the cor- 
responding relic abundance is typically smaller than the estimated non-baryonic matter 
content of the Universe. This is first due to the direct annihilation rate of higgsino-like 
or wino-like LSPs, which is by far larger than that of a bino-like LSP (see e.g. the left 
panel of Fig ^). In fact, not only the couplings involving, for instance, a wino are larger 
than those involving a bino, since g2 > Qi, but also the number of final states in which 
winos or higgsinos can annihilate is larger. As an example, winos can annihilate into a cou- 
ple of weak bosons without any intermediate scalar superpartners, whereas binos cannot. 
Furthermore, the mass spectrum of a wino-like LSP always yields an approximate mass 
degeneracy with the lightest chargino, as does a higgsino-like LSP, which is, in its turn, 
also very close in mass to the next-to-lightest neutralino. This means that coannihilations 
between the LSP and the lightest chargino, as well as, in the case of the higgsino, with the 
neutralino X2, further reduce the relic abundance in the whole parameter space. This mass 
degeneracy with the chargino also bears a rather strong lower bound on the neutralino 
mass through direct accelerator searches. As a result, in the absence of relic density en- 
hancement processes^, wino or higgsino LSPs are not adequate to be the main dark matter 
component. 

In the present framework of quintessential cosmologies, the enhancement resulting from 
the shift towards larger values of the freeze-out temperature of cold relics may naturally 
render wino or higgsino LSPs attractive dark matter candidates. In order to elucidate 
this point, we will discuss two particular cases: mSUGRA with non-universal gaugino 

Among these, we mention non-thermal neutraUno production through direct or indirect decays of 
gravitinos or moduh fields p^]. 
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masses generated by a gauge kinetic function belonging to the 200 representation of the 
symmetric product of two SU{5) adjoints, where the LSP is mostly higgsino-like, and the 
minimal Anomaly Mediated SUSY Breaking (mAMSB), where the LSP is typically almost 
completely a wino. We show that (1) the relic abundance naturally falls into the preferred 
cosmological range for a suitable choice of the quintessential parameters and (2) that the 
needed fine tuning in the cosmological parameters is reasonably low, and comparable with 
that of the supersymmetric parameters. 

4.1 SU{5)-like gaugino non-universality 

In the context of supergravity, non-vanishing gaugino masses are generated through the 
SUSY breaking vev of the auxiliary component of the lowest order non-renormalizable 
term in the gauge kinetic function. The SUSY breaking vevs, in order to preserve gauge 
invariance, must lie in a representation of the symmetric product of two adjoints of the 
underlying unified gauge group. In the case of SU{5) GUTs, one has 

(24 X 24)3y^„^ = 1 © 24 © 75 © 200. (4.1) 

Universal gaugino masses are generated only if the SUSY breaking field lies in the singlet 
representation 1. In all other cases, particular ratios between gaugino masses will hold at 



the GUT scale |25|, We report in Tab. y the resulting relative gaugino masses p,8 \. 

Two of the four representations yield a higgsino-like LSP, namely the 75 and the 
200. We choose here to focus on the case of the 200, because in the case of the 75 
successful electroweak symmetry breaking (EWSB) forces the parameter space to rather 
narrow regions, and to low values of tan/3 p^ . 

We study the requirements on the cosmological sector needed to achieve a neutralino 
relic density compatible with the current estimates of the dark matter content of the 
Universe. In particular, we consider a sample case in which we fix tan f3 = 45, Aq = 0, 
/X > and mo = 500 GeV, and let M1/2 vary within the range allowed by EWSB. Changing 
the value of tjiq does not lead to remarkable changes in the results: we only find that at lower 
values TJiQ < 500 GeV and at low Mi ^2 the j4- funnel condition ~ niji/2 approximately 
holds, while at larger values of mg the allowed M1/2 range is more and more reduced by 
the requirement of successful EWSB. In order to have a wider range of neutralino masses, 
we resort to a large value of tan /3. 



Rep. 
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24 


-1/2 


-3/2 


1 


75 


-5 


3 


1 


200 


10 


2 


1 



Table 2: Relative values of the gaugino masses at the GUT scale for different representations 
allowed in SU{5). 



-17- 




150 200 250 300 350 400 450 500 150 200 250 300 350 400 450 500 

ni/GeV) mj(GeV) 

Figure 9: {Left ): The (m^, Hi) plane for the case of the representatfon 200, with SUSY param- 
eters tan/3 — 45, Aq — 0, n > 0, niQ — 500 GeV and fixing A = 3.5. The yellow shaded region has 
fi^/i^ > 0.13, while in the green-shaded strip 0.09 < Q^h'^ < 0.13. {Right ): The (m^, plane 
at the same values of the SUSY parameters. In both panels, the blue dotted lines correspond to 
points at fixed nh'^ = 0.005, 0.01, 0.025, 0.05, 0.09, 0.113, 0.13, while the black dashed lines to 
points at fixed enhancements Afi — 1, 5, 10, 25, 50. 

We report our results in Fig where rather than M1/2 we consider the neutralino mass 
as the last SUSY parameter which fully fixes our model setup: the mass range displayed 
is dictated by the fulfillment of phenomenological bounds and of successful radiative EWSB. 
The Quintessence contribution to the expansion rate of the Universe is defined by fixing 
the parameter A = 3.5 in the exponential potential eq. ( p.5D and by varying the initial 
condition on Hi. In the left panel we show the parameter Hi, while in the right panel we 
plot a more model independent result by trading Hi for the Quintessence-to-radiation ratio 
at neutralino freeze-out The region shaded in yellow bears an excessive relic density 
enhancement, giving rise to > 0.13, while the green shaded strip has 0.09 < ^^^^^ ^ 

0.13, and thus approximately reproduces the 2-cr range for f^cDM^^- In the figures we 
also plot some isolevel curves both for O/i^ (blue dotted lines) and for AO, (black dashed 
lines). At large Hi, i.e. when Quintessence is dominant only at very early times, the 
iso-il curves tend to become vertical lines, in the region where AQ <ti 1. As expected, we 
find that the quintessential enhancement, quantified by AQ, is, to a good approximation, 
fixed by as emerging from the (m^,^^) plane, where the isolevel curves for AO, are 
practically horizontal. The peculiar behavior of the iso-fi curves as functions of is 
due to the fact that, raising m^, the higgsino fraction of the lightest neutralino grows, 
thus partly compensating (at < 250 GeV) the effect of a larger mass on the relic 
density. When the neutralino starts to be largely higgsino-dominated (m^ > 250 GeV), 
monotonously grows with the neutralino mass. We highlight that at the left sides of 
the figures, i.e. at low m^, the required enhancement is AO ~ 50, which means that the 
neutralino relic density in the absence of Quintessence is two orders of magnitude below 
QcDM- Actually, in the absence of Quintessence, the neutralino relic density, at the smallest 
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Figure 10: {Left ): The {m^, Hi) plane for the case of the anomaly mediated SUSY breaking, 
at tan/3 = 50, /i > and mo — 1 TeV. The yellow shaded region has fl^^h^ > 0.13, while in 
the green-shaded strip 0.09 < fi^/i^ < 0.13. {Right ): The (m^, plane at the same values 
of the SUSY parameters. In both panels, the blue dotted lines correspond to points at fixed 
nh'^ = 0.0005, 0.001, 0.005, 0.01, 0.025, 0.05, 0.09, 0.113, 0.13, while the black dashed lines to 
points at fixed enhancements Ail = 0.1, 1, 5, 10, 25, 50. 

value for is found to be Qh'^ ~ 0.002. We checked that the 200 model yields, at any 
tan/3, a neutralino relic density lying at least a factor 2 below Ocdm) and in the specific 
case tan/3 = 45.0, f^^^^ ^ 0.03: therefore, in this context, a relic density enhancement 
mechanism is mandatory. 

4.2 mAMSB and wino LSP 

In the minimal anomaly mediated supersymmetry breaking (niAMSB) model, the soft 
SUSY breaking contributions arising from the super- Weyl anomaly are supposed to be the 
dominant ones [^0| |3l]] . As a consequence, the gaugino masses are entirely determined by 
the corresponding /3 functions and by the gravitino mass rn-3/2- The problem of negative 
scalar masses is cured by adding a universal contribution m-Q. The model is then defined 
by a set of four parameters, namely 

mo, ""1.3/2, tan /? and sign^u. (4-2) 

Since at the weak scale M1/M2 ~ 3.2, the LSP of mAMSB models is dominantly a wino, 
with a very large typical purity. The lightest chargino is thus quasi degenerate with the 
LSP, and the total cross section of neutralino pair annihilation is very large The 
pattern for the quantity W^s is therefore rather similar to that of the higgsino LSP case, 
a part from the absence of the coannihilations with the next-to-lightest neutralino, and 
therefore of the second bump visible in the left panel of Fig ^. As far as the quintessential 
enhancement is concerned, we find, as expected, that the enhancement pattern for wino-like 
neutralinos is very close to that in the higgsino case. 

We study in Fig ^ the quintessential relic density enhancement needed to obtain 
cosmologically viable neutralino relic densities in the wino LSP mAMSB scenario. We 
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Figure 11: The fine timing witli respect to tlie cosmological parameters A (red solid lines) and 
log[-ffi] (dotted blue lines), and to the supersymmetric parameter tan/3 (dashed magenta lines). 
With dot-dashed green lines, we also show the fine tuning due to Mi ^2 i^wf^i) ^'^d to to 3 (left), 
respectively (right) in the case of higgsino dark matter in mSUGRA with non- universal gaugino 
masses (rep. 200) and {lejfj in the niAMSB scenario. The amount of fine timing is computed along 
the curves at fixed flh^ = 0.113, as a function of the neutralino mass m^. 



focus on a sample case at mo = 1 TeV, tan/3 = 5, fi > and ranging from 32 to 192 
TeV. We plot in the figures, as in the previous higgsino case, the neutralino mass range 
allowed by accelerator searches (lower bound) and by successful EWSB (upper bound). 
Again we shade the cosmologically preferred region in green, while the yellow shaded parts 
yield an overproduction of neutralinos. We plot the {rriy., Hi) plane in the left part of 
the figure, while {m^, 1^^) in the right panel. Also shown are the curves corresponding to 
particular values of the resulting relic density (blue dotted lines) and of the quintessential 
enhancement (black dashed lines). As emerging from the comparison between Fig ^ and 



Fig |10|, in the present mAMSB case, the needed enhancement factors are larger than in the 
higgsino LSP case, since wino pair annihilations have a larger cross section with respect 
to the higgsinos. In mAMSB scenarios, the wino-purity does not critically depend on the 
neutralino mass, thus the iso-Jl curves are remarkably smooth. The relic density, in its turn, 
is a growing function of m^, and therefore the quintessential enhancement needed to drive 
the resulting neutralino relic abundance to the cosmologically preferred range decreases 
with increasing neutralino masses. We remark that in the present case the neutralino 
relic density in the absence of Quintessence is always very low, being ~ 0.0003 in 



correspondence to the lowest neutralino masses. In fact, in the case shown in Fig 10, the 
needed enhancement factor can be as large as Ail ~ 300. 

4.3 Fine Tuning 

The issue of fine tuning in supersymmetric models of dark matter [^] is often quantified 
through the logarithmic sensitivity of the neutralino relic density with respect to the various 
parameters appearing in a given model. The amount of fine-tuning for a given phenomeno- 
logical model is often regarded as a measure of its degree of naturalness. Moreover, the 
experimental accuracy needed to fix a given parameter is, typically, inversely proportional 
to the sensitivity of the measured quantity on the parameter itself. It is then of crucial 
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importance to understand whether quintessentially enhanced models of dark matter are 
"natural" or excessively fine-tuned. We therefore concentrate, in this section, on the fine 
tuning in the cosmological sector of a supersymmetric model embedded in a quintessential 
cosmology. 

To estimate the fine tuning in the parameters defining the cosmological model, i.e. Hi 
and A, following we compute the quantity 



We show our results in Fig. 11. The amount of fine-tuning is quantified in the two cases 
discussed in the preceding sections 4.2 and in correspondence to the lines at Qh'^ = 



0.113, i.e. the central value of the dark matter density as indicated by the latest evaluations 
|T|. The red solid lines refer to the case aj = A, while the blue dotted lines to = 
log Hi. For comparison, we also include the amount of fine-tuning in the supersymmetric 
parameters tan P and M1/2 (^-3/2) foi' the 200-higgsino like neutralino (resp. mAMSB). We 
notice, as expected, that the most fine-tuned parameter is Hi, which, in the present context, 
sets the amount of Quintessence in the early Universe. The dependence on A is, on the 
other hand, found to be not critical, yielding a moderate fine-tuning in both the higgsino 
and wino dark matter cases. As compared with the supersymmetric parameters taken into 
account, we find that the fine-tuning in Hi is of the same order of magnitude of that of 
tan P in the higgsino-like neutralino case, while in the context of mAMSB the dependence 
on tan /3 is largely suppressed: this is traced back to the fact that in the higgsino case, both 
Ml/2 and tan/? set the higgsino fraction in the lightest neutralino, and therefore, due to 
the large difference in the pair annihilation of a bino and of a higgsino, they play a crucial 
role in determining Qh"^. The sensitivity of the neutralino relic density to the quintessential 
parameters is thus found to be of a certain relevance only for the parameter Hi, though the 
corresponding fine-tuning may be of the same order as that needed in the supersymmetric 
parameters in the sample models considered here or even in the standard mSUGRA setup: 
in particular even larger amounts of fine tuning are in order in the focus-point, funnel or 



sfermion coannihilation regions, see, e.g., |34, 



5. Conclusions 

We analyzed the quintessential enhancement of the neutralino relic abundance within "re- 
alistic" cosmological models. We chose to resort to models of Quintessence exhibiting 
tracking solutions stemming from an exponential potential, and characterized by only two 
parameters: the exponent appearing in the potential and the Hubble rate at the end of 
inflation. In this case, we numerically computed the evolution for the expansion rate pa- 
rameter and interfaced the result with a fully-numerical computation of the neutralino 
relic abundance, with an overall estimated accuracy of the order of 1% or better. We 
defined a suitable parameter to quantify the enhancement in the relic abundance induced 
by Quintessence and found the limit above which a Quintessence component becomes rel- 
evant: non-negligible effects are expected, at the 1% or higher, whenever the ratio in the 
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Quintessence to radiation energy densities at neutralino decoupling is larger than about 
0.1. Furthermore we showed that, within our viable cosmological models, the enhance- 
ment in the neutralino relic density can be as large as six orders of magnitude. Then, we 
systematically analyzed the dependence of the enhancement on the features of the lightest 
neutralino as a dark matter candidate (such as its mass and composition or any particular 
mechanism, e.g. coannihilations or resonance effects, playing a role in the neutralino relic 
abundance), concluding that the largest effects take place at large masses and in the coan- 
nihilation regions. We also showed that models with a dominantly wino- or higgsino-like 
lightest neutralino (as it the case in mAMSB or in SU (5)-inspircd gaugino non-universality 
with the 75 or the 200 representations) may well be compatible with the required dark 
matter content of the Universe if a quintessential enhancement is present. We also showed 
that, in the proposed models, the fine-tuning in the cosmological parameters is always 
rather low. 
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